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ABSTRACT 

The intensity mapping of the [C ii] 157.7 fine-structure emission line represents an 
ideal experiment to probe star formation activity in galaxies, especially in those that 
are too faint to be individually detected. Here, we investigate the feasibility of such an 
experiment for z > b galaxies. We construct the Ten ~ -^h relation from observations 
and simulations, then generate mock [C ii] intensity maps by applying this relation to 
halo catalogs built from large scale N-body simulations. Maps of the extragalactic far- 
infrared (FIR) continuum, referred to as “foreground”, and CO rotational transition 
lines and [C i] fine-structure lines referred to as “contamination”, are produced as 
well. We find that, at 316 GHz (corresponding to zcii = 5), the mean intensities 
of the extragalactic FIR continuum, [C ii] signal, all CO lines from J = 1 to 13 
and two [C i] lines are ~ 3 x 10® Jy sr“®, ~ 1200 Jy sr“®, ~ 800 Jy sr“^ and 
100 Jy sr“^, respectively. We discuss a method that allows us to subtract the FIR 
continuum foreground by removing a spectrally smooth component from each line 
of sight, and to suppress the CO/[C i] contamination by discarding pixels that are 
bright in contamination emission. The z > b [C ii] signal comes mainly from halos in 
the mass range 10^®“^^ as this mass range is narrow, intensity mapping is an ideal 
experiment to investigate these early galaxies. In principle such signal is accessible to 
a ground-based telescope with a 6 m aperture, 150 K system temperature, a 128 x 128 
pixels FIR camera in 5000 hr total integration time, however it is difficult to perform 
such an experiment by using currently available telescopes. 

Key words: cosmology: diffuse radiation-dark ages; reionization, first stars - radio 
lines: galaxies - galaxies: high-redshift 


1 INTRODUCTION 


The emission line intensity mapping is a technique to ac¬ 
cess high-z galaxies below the detection limit without los- 


ing redshift information, as 
(l201(tl: Visbal et all ll201lh: 

Droposed by e.e. IVisbal & 
Gone et al. 1I2OIII. 20ll 

c Loeb 

20121: 

iLidz et aLl 1I2OIII): ISilva et al.l (20131): iPullen et al.l 1 

20141. 


Optimistically, it only collects radiation from galaxies in 
a selected redshift range, as the spurious flux due to fore¬ 
grounds, contaminating radiation and noise can be in prin¬ 
ciple removed or suppressed. Compared with galaxy surveys 
that aim at resolving faint spots in a limited field of view 
(FOV), the advantages of intensity mapping rely on the fact 
that, if the galaxy luminosity function has a sufficiently steep 
faint end, the obse rved radiation is actually dominated by 
unresolved sources (lUzgil et al.ll20lj l. Even if this is not the 
case, intensity mapping can still be used to study unresolved 
galaxies once resolved sources are removed (masked). Inter¬ 


estingly, an intensity mapping experiment could be carried 
out with a modest aperture but large FOV telescope. 


The [C ii] 157.7 /rm fine-structure line arising from 
3/2—>-^Pi/2 transition is the brightest amongst all 
metal lines emitted by the interstellar medium (ISM) of 
star-forming gala xies. It is associat e d to the star forma¬ 
tion in galaxies (|Boselli_gt-aL 200^ de^JjOozeeUal I I2OIII : 


iDe Looze et al.ll2014l : lHerrera-Camus et al.ll2014| 'l and plays 

a key role in the energy balance of galaxies, as it provides 
one of the most efficient cooling processes for the neutral 
ISM. With respect to the Lya line, the [C 11] line has the 
advantage of being unaffected by dust attenuation and neu¬ 
tral hydrogen absorption. 

In the local Universe, [C ii] line has been success¬ 
fully detected even in galaxies with amazin gly low star for- 
mation rates (SFR) of ~ 0.001 MQyr~^ JPe Looze et al.l 
|2014| . These authors have also derived the relation be¬ 
tween the [C ii] line luminosity, Lcii, and the SFR of local 
galaxy samples (IPe Looze et al.ll2014l f. Surprisingly, given 
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the rather complicated physics behind the [C ii] emission 
Leu scales rather tightly with SFR. However, at high red- 
shift {z ^ 4), the [C i i] line has been detec t ed so far only in 
quasar host galaxies |Maiolino et al. 200^ Vmiemmis_et_^ 


2012 


2013 


Gallerani et al.ll2012l : Wang et al. 20131 : Willott et al 

Cicone et alJ 12011 7 or ultra-luminous infrared galax¬ 


ies (ULIRGs, with Lib. > 10^^ L© where Lib is the in- 
band luminosity a t 8 — 1000 /rml characterized by SFR 
~ 10^-3 MQyT~^ lICox et al.ll201ll : iDe Breuck et al1l201ll . 
I2OI4I I. For typical normal star-forming galaxies (SFR ~ 
10 M©yr~^) , [C I II emission has not yet been detected 
llOuchi et al'l 1^131: lOta et al.l l2014l : ISchaerer et alj l2014l : 
iGonzalez-Lopez et al. 2014l ~l. This might indicate that most 
of carbon in these galaxies is at higher ionization state 
and/or their ISM is characterized by a very low level of 
metal enrichment. By applying the Lcii—SFR relation de¬ 
rived from local galaxies samples to high redshift Lyo emit¬ 
ters it is possible to compute the expected [C ii] flux from 
these galaxies. The fact that their [C ii] line remains unde¬ 
tected even with ALMA provides useful constraints on their 
internal rad iation field, molecular content, gas density, and 
meta llicity dVallini et al.l[2015l . l2013l : iGonzalez- Lopez et al.l 

I2OI4IL 


As probing tools, intensity mapping experiments are af¬ 
fected by the presence of foreground radiation, including 
that represented by the galaxy continuum redshifted into 
the observed band. Unfortunately, it is almost often the case 
that the foreground intensity largely exceeds that of the sig¬ 
nal. The typical [G ii] lin e luminosity is 0.1% — 1% of the 
Lib (iMaiolino et al .I2OO9IL This implies that even if only one 
percent of the IR luminosity is redshifted into the observed 
band, the continuum emission overcomes the [G ii] line. 

In addition to the far-infrared (FIR) continuum fore¬ 
ground, there are other emission lines emitted from a range 
of redshifts that fall at the same frequency of the [C ii] 
signal; they act as contaminants. For example, the [O i] 
line with wavelength 145 pm, the two [N ii] lines (A = 
122,205 pm) and two [C i] lines (A = 610,371 pm), and 
a handful of GO rotational transition lines in the range 
200-2610 pm. Among these, the GO rotational transition 
lines are the most relevant here. For example, since the 
GO(4-3) line has a wavelength 651 pm, if emitted from 
2 = 0.45 galaxies, it contaminates the [G ii] emission from 
z = 5 galaxies. The emission efiiciencjQ of the CO (4-3) 
line from star-formin g galaxies is ~ 2% of the [C ii] line 
dVisbal fc LoebI I2OICII L However, the luminosity distance 
from 2 = 0 to 0.45 is only ~ 5% of that to 2 = 5. As 
the flux is inversely proportional to the square of the lumi¬ 
nosity distance, whereas the proper distance interval that 
corresponds to the same bandwidth is oc (1 -I- 2)“®^^, the 
CO flux can be more than ten times higher than the [C ii] 
one, even ignoring the cosmological evolution of the star 
formation rate density. Thus, CO contamination, as well as 
the continuum foreground, cannot be ignored and must be 
considered thoroughly. 

Alth ough the (C ii] signal itself can be computed ana¬ 
lytically dCong et al.ll201^IUzgil et al.ll2014h . a reliable in¬ 
vestigation of the influence of foreground/contamination is 


only possible based on mock maps that carefully mimic ob¬ 
servations as close as possible. This is the prime motivation 
of this paper. Using halo catalogs recovered from large scale 
N-body simulations, we produce mock maps including (a) 
[C ii] signal, (b) FIR continuum foreground, (c) CO and 
[C i] contamination lines, and (d) instrumental noise. We 
then test our foreground/contamination removal scheme on 
these maps to demonstrate the successful recovery of the 
original [C 11] signal. 

The layout of this paper is as follows. In Sec. O we 
describe the model used to compute the [C ii] from high-2 
galaxies and the necessary steps to generate mock maps. 
We show our forecasts for the [C ii] signal, extragalactic 
FIR continuum and the contamination and perform fore¬ 
ground/contamination removal experiments on mocks to re¬ 
cover the original [C ii] signal. The results are presented in 
Sec. [3] The conclusions and discussion are found in Sec. ID 


2 METHODS 

2.1 [C ii] emission from early galaxies 

IVallini et al.l ll2013h and IVallini et al.] (l2015l ) (hereafter V15) 
have combined high-2 galaxy numerical simulations with 
sub-grid models of the ISM to compute the expected [C ii] 
luminosity (Lcii) arising from diffuse neutral gas and pho- 
todissociation regions (PDRs). The resulting trend of Lcii 
with SFR and metallicity (Z) is consi stent with observation s 
of local metal-poor dwarf galaxies 1^ Looze et al.l l2014l L 
For the range [0.1,100] MQyT~^ and [0.05,1.0] Zq which 
should encompass most of the sources contributing to the 
total ]C ii] emission at high redshift, V15 results are well 
reproduced by the following fitting formula 

log(Lcii) = 7.0 + 1.2 X log(SFR) -f 0.021 x log(Z) 

+ 0.012 X log(SFR)log(Z) - 0.74 x log^(Z), (1) 

where Lcii, SFR and Z are in units of L©, M©yr“^ and Zq 
respectively. 

The next step is to compute the Lcii — Mh relation, 
where Mh is the halo mass. To this aim we need to know 
the SFR—Mh and Z — Mh relations. We obtain them from 
the procedure described in below paragraph^. 

Since the UV luminosity (Luv) of a galaxy scales with 
its SFR fe.g. lKennicutt|[l998l L we adopt the observed UV 
luminosity functions (LFs) to derive the SFR—Mh relation. 
The measured UV LF is well d escribed by a Schechter pa¬ 
rameterization (ISchechteilll976tl : 

^ = 0.41n(10)^*.-“e-, (2) 

where x = lo'^ '^^^uv-^uv) ^ with Muv the dust-attenuated 
absolute magnitude. For the rest-frame UV luminosity at 
1600 A, the redshift-dependent parameter s 
that fit observations between 2^4 — 8 are (|Bouwens et al.l 


^ As a caveat, we note th at there is no clear consensus in the 
literature on this value, see lBrevsse et al] ll2014lV 


^ Note that, compared to galaxies , the inter^alacti c medium 
emits a negligible [C ii] signal (e.g. iGong et al.l 12013) at z ^2^ 
and therefore is not considered in this work. 
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l2014bl) 

Mw = -20.96 + 0.01(z - 6) 

0* = 0.46 X (3) 

a = -1.87- 0 . 10 ( 2 -6). 

The intrinsic absolute magnitude is Myv = Afuv — ^leoo, 
where A iann = 4.43 + 1.99/3 (^0) is the dust attenuation at 
1600 A dMeurer et al.l[l999ll and /3 is the measured spectral 
slo pe (/a oc A^). Normally p depends on Muv and is fitted 
by (iBouwens et al.|[2014af l 

/3 = /3-i9.5 + ^^(Muv + 19.5). (4) 

From Fig. 2 of iBouwens et af] (l2014al l we find the following 
redshift-dependent fit, valid for A P, z P, 7 


P-1Q.5 

dp 

dMjj-v 


-1.97- 0.06(2 - 6) 


= -0.18-0.03(2-6). 


(5) 


The intrinsic UV LF is then connected to the measured UV 
LF via 


dm! 

dA^;^ 


(M(jv, 




dn 

dM\jv 


{Mvv, z). 


( 6 ) 


Assuming that the intrinsic Tyv monotonically in¬ 
creases with Mh and that all halos host some star formation 
activity, we obtain the Tuv — Afh relation from 



dn' 

d^ 


dMuv = 



dn 

dM^ 


dMh, 


(7) 


where dn/dMh is the h alo mass function (ISheth fc Tormeiil 
1 19991: ISheth et al. ][20 o 3). This “abundance matching” tech¬ 
nique will be used also in Sec. m\ to derive the relation 
between the IR luminosity and halo mass. We then derive 
the SFR from Lyv In principle Tyv depends not only on 
the SFR, but also on metallicity and stellar age. However, we 
note that the UV luminosity is insensitive to the metallicity 
and the stellar age, unless stars are very young ( ,$ 10 Myr). 
So we can safely assume that Lyv scales with SFR as 


L{jy = Ivv X SFR. (8) 

We compute Ijjy from Starburstgj^ (iLeitherer et al. I l 19991: 

I Vazquez fc Leithererl[2005l : ILeitherer et al.l 2010l 'l bv" assum¬ 
ing a metallicity 0.1 ^0, stellar age 10% of Hubble time, 
and a Salpeter IMF between 0.1 — 100 Mq. We choose 
the “continuous star formation” mode. At 1600 A, lyy = 
(8.9, 8.6, 8.3) X 10^ ^ erg s~^Hzy (M pi/yr)"^ at 2 = (5,6,7); 
which is similar to iKennicuttI lll998li . 

We plot the SFR derived with this procedure at various 
redshifts as a function of Mh in Fig. [T] As a co mparison we 
also p lot the SFR—Mh relation at 2 ~ 5 found bv ISilva et al.l 
(l2014h who fit the relation using semi-analytical models of 
galaxy formation. As can be seen by inspecting this figure, 
the specific SFR, i.e. the SFR per unit mass, starts to drop 
at a turnover mass ~ 10^^ Mq. This is consistent with semi- 
analytical model predictions. 

The final ingredient of Eq. © is Z. As the metallicity 
of high -2 galaxies is very poorly constrained at present, we 


® http://www.stsci.edu/science/starburst99/docs/default.htm 


derive it by combining the Lyv — Mh relation with the “fun¬ 
damental metallicity relation” (FMR) that relates Z to the 
stellar Mass (M*) and SFR. The FMR inferred from low -2 
galaxy observations (iMannucci et al.ll20ldf is given by the 
following equation: 


log(Z) = 0.21 -h 0.37 log(Mio) - 0.14 log(SFR) 

- 0.191og^(Mio) - 0.054 log^(SFR) (9) 

-f 0.121og(Mio)log(SFR), 

where Mio = M*/10^°, M* and Z are expressed in Solar 
units. No redshift evo lution is found at least up to 2 = 2.5 
llMannucci et al.ll2010lf : therefore we apply it also to high -2 
galaxies, with the caveat that deviations might appear for 
high -2 galaxies. The stellar mass, M*, is linked to the UV 
absolute magnitude via the mass-to-ligh t ratio. From the 
latest measurements (iDuncan et al ][2oii), 

log(M*) = log(Mr) + + 19.5), (10) 

where log(M*) = (9.00,8.84,8.63) and dM*/dMuv = 
(—0.46, —0.54, —0.45) in the redshift ranges 4.5 ^ 2 < 5.5, 
5.5 ^ 2 < 6.5, and 6.5 ^ 2 < 7.5, respectively, for the model 
without nebular line contribution. We use the observed UV 
magnitude Muv (i.e., without dust correction), and we ig¬ 
nore the negligible difference between the UV luminosity at 
1500 A and 1600 A. Fig.[2]shows the derived Z — Mh relation 
at 2 = 5, 6 and 7. 

We finally compute the [C ii] luminosity of halos with 
mass Mh by substituting the SFR and Z derived above 
into Eq. ©• Although different star formation histories may 
cause a scatter in the [C ii] luminosities for halos of a given 
mass, we neglect this effect because it only results in noise 
as long as the luminosity dispersion is independent of posi¬ 
tion on large scales. The Lcii — Mh relation derived with the 
above procedure is shown in Fig. [3] by solid, short dashed 
and long dashed lines for 2 = 5, 6 and 7 respectively. 

When deriving the above Lcii — Mh relation, the prop¬ 
erties of faint galaxies are extrapolated from the observed 
bright galaxies. We check the validity of this relation by 
comparing the result obtained through this semi-empirical 
method with the numer ical simulations of co smic metal 
enrichment presented by IPallottini et al.l ll2014h (hereafter 
P14). P14 have used an hydrodynamical simulation to fol¬ 
low the star formation and the Pop Ill-Pop II transition 
{Z > Zq). Hereafter, we only consider the Pop II star 
formation mode. 

If a halo has formed stars, it contains a number of stellar 
particles whose birth date and metallicity is recorded. For a 
selected halo, the mean stellar age is 


t 


age — 




( 11 ) 


where Ati is the birth date of the i-th stellar particle, m-t^i 
is its mass and the sum is extended over all Pop II stellar 
particles in that halo. The mean metallicity of the halo is 


Z = 




( 12 ) 


where Zi is the metallicity of the i-th Pop II stellar particle. 
By dividing the total Pop II stellar mass by the mean age. 
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Figure 1. The SFR derived from UV LFs as a function of halo 
mass at 2 : = 5, 6 and 7 respectively. The SFR—Mh relation in 
ISilva et al.l ll2014l') at 2 : 5 and the SFR of each halo containing 

Pop II stars in P14 simulation at 2 : 5 are also shown. 



Figure 2. Metallicity, derived from mass-to-light ratio and the 
FMR, as a function of halo mass at 2 : = 5, 6 and 7. The metallicity 
of each halo containing Pop II stars in the P14 simulation at 2 5 

are also shown. 


we obtain the mean SFR of the halo, 

SFR= (13) 

^age 

The SFR vs. Mh and Z vs. Mh at z ~ 5 in P14 simulation 
are also plotted in Fig. [T]and Fig. [2] respectively. 

Using Eq. m we calculate the [C ii] luminosity of halos 
from their SFR and Z and group halos into several mass 
bins at each simulation output. Some halos only have Pop 
III stars, or are too small to host any star formation. There¬ 
fore in each mass bin only a fraction feu of halos exhibit 
[C ii] emission. This fraction tends to one as the halo mass 
increases. We denote the mean log of [C ii] luminosity for 
halos that exhibit [C u] emission by (log(Lcii)), and use 

log(Tcii) = (log(Lcii)) -f log(/cii) (14) 

as the mean [C ii] luminosity of all halos with mass Mh; this 
quantity is plotted in Fig. [S] 

2.2 Far-infrared continuum foreground 



M, [Ma] 


In this subsection, we model the extragalactic foreground 
due to the FIR continuum from galaxies at different red- 
shifts. The Milky Way FIR and CMB radiation are assumed 
to be removed straightforwardly and hence are not consid¬ 
ered in this work. 

The FIR luminosity function of galaxies, includ¬ 
ing spiral galaxies, starburst galaxies, star-forming galax¬ 
ies contai ning AGNs and sometimes AGNs, i s stud ¬ 
ied in e.g. Magnelli et al.l ll2009l'l: Gruppi oni et al.l (l2013ll : 
iMagnelli et al.l ( 2013lf . In iGruppioni et al.r i 2013ll . the LF 


can be written as 


exp 


1 1 2 


i-b 


Tir \ 

LmJ\ ’ 


(15) 


where Lir is the infrared luminosity between 8 — 


Figure 3. The mean Lqu ^ function of halo mass, Mh, at red- 
shifts ~5, 6 and 7 respectively. The lines are the relation derived 
using observations while the points are for the P14 simulation 
(the errorbars are the standard deviation of {log(Lcii)) for each 
mass bin). 


1000 fim. We use redshift evolution formulae of parameters 
a,a, (jCruppioni et al.ll2013h : = 5.7 x 10“^(1 + 

^)-o.57 ^ ^ ^ ^ 10"^(1 + for 

2 > 1.1; = 7.68 X 10^(1 + for 2 ^ 1.85 and 

L^r = 5.80 X 10^'^(l + z)^-®^ for 2 > 1.85; a - 1.15, a = 0.52 
for z ^ 0.3, and a = 1.2, a — 0.5 otherwise. The above LFs 
are constructed from galaxy samples at z ^ 4.2, therefore 
including the large majority of the sources contributing to 
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Figure 4. FIR continuum luminosity derived from Eq. Cii as a 
function halo mass at redshift 2 = 0.5 and 2. 


we assume a Gaussian distribution of the form 


p(Lco| Leo) = 


V27r( 


-exp 


X 


(17) 


where x = log(Z/co) — log(ico) and a = s/a is the variance; 
s is the s catter of the data around the log(I/iR) — log(ico) 
htting in iGreve et alj (l2014ll . including the intrinsic disper¬ 
sion and statistical errors. 

We further consider the contamination from two [C i] 
fine-structure lines: (i) [C i(l-O)], corresponding to the 
®Pi—^-^Po transition, at 492 GHz, and (ii) and [G i(2-l)], 
corresponding to the ®P 2 —^®Pi transition, with frequency 
809 GHz. Seve ral authors have reported [C i] observa¬ 
tional data (e.g. Gerin fc Phillips! 200^ Israel fc Baasllioo^ : 
IWalter et alJ 1201 ll : Alaghband-Zadeh et al.l 20131 '). finding 
relations between the [C i] and CO or IR lum inosities. Moti¬ 
vated by observations. [Popping et al.l ll2014h have calculated 
the expected [G i]/I/ir ratios for 0 < 2 : < 2 galaxies, by 
combining a semi-analytical galaxy formation model with 
radiative-transfer and line-tracing calculations. We adopt 
the outcome of these theoretical calculations and we add 
0.25 dex scatters to the mea n ratio, namely the m aximum 
of the deviations reported by [Popping et all ([2014[ i. 


the FIR continuum (and GO, which is associated to the FIR 
continuum, see next subsection). 

We use again th e abundance matching technique 
llBethermin et aL[[2012[ l to construct the Lir — Mh relation. 
We suppose that the contribution of subhalos to the IR lu¬ 
minosity function is small and we ignore them. By equating 
the number density of galaxies with IR luminosity above Lir 
and the number density of halos above Mh, 


/ '1 ?(Lir, «)dLiR = / 

J LiYi -'Afh 


dn 

dAP) 


dMi,, 


(16) 


the Lir — Mh relation is derived. 

We plot the the IR luminosity - halo mass relation at 
redshift 0.5 and 2 in Fig. [I] For the same reasons given in 
Sec. 12.11 we do not consider the IR luminosity dispersion 
among halos with the same mass Mh. 


2.3 CO and [C i] emission lines contamination 

The CO rotational transition lines from low-redshift galax¬ 
ies are by far the most important contaminants for the 2 > 5 
C iisignal. Although some stud ies aimed at mea suring CO 
luminosity functions exist (e.g. [Keres et al. | [2003l and refer¬ 
ences therein), they lack data either for higher rotational 
transition numbers, J, or at high redshift. On the other 
hand, the CO line luminosity is found to be closely related to 
the IR luminosi ty, since both line s are good star forma tion 
activity tracers l|Bavet et al.[[200^ : [Popping et al.l[2014[ 'l. 

The CO line luminosity is derived by using the 
log(LiR) = Q log(L co) -|- P relations presented in Tab. 3 of 
[Greve et al.l ([20141 1 for lines with J-ladders from 1 to 13. 
These relations are fitted from samples of local (2 < 0.1) 
(U)LIRGs and high -2 (2 > 1) dusty star-forming galaxies 
(DSFGs). Since CO lines are considered as contaminants to 
be removed when recovering the [C ii] signal, the scatter 
around the fitting relation should be modeled. To this aim 


2.4 Instrumental noise 


We have to account for instrumental noise, in order to have 
predictions that can be fairly compared with observation. 
The noise level of a radio telescope is given by the standard 
expression 


ctn 


2kBTaya 

A\/Avot ’ 


(18) 


where ks is the Boltzmann constant, Taya is the system tem¬ 
perature, A is the area of the antenna, t is the integration 
time per FOV. If the camera has Apix pixels and the observa¬ 
tion is performed at wavelength Ao, normally the instrument 
is designed by such a way that Hfov ~ Ap/AApix holds. To 
cover a sky region of solid angle Hmap, the total integration 
time is fobs = t x Hmap/HFOV- We model the instrumental 
noise as a zero mean Gaussian random variable without spa¬ 
tial and frequency correlation, and we add such fluctuations 
to the mock maps. 


2.5 Mock maps 

The light cone for which we produce the intensity maps is 
built from the halo catalogs of the BolshoiP simulatiorQ 
(Bolshoi simula tion with Planck c osmology, see Bolshoi 
simulation paper lxiypin et al.l[201l| l. In the simulation, the 
smallest halos resolved are « 5 x 10® Mq, well below the 
mass of halos that are expected to host the bulk of [G ii] 
emission (see Sec.s 1 ^ . 

The box of this simulation is L = 250 h“^cMpc on 
a side, corresponding to 2.4 degree when located at 2 = 
7. When making light cones from a simulation with a box 
length smaller than the cone depth, a standard protocol is to 
replicate the same halo catalog along the radial direction at 
the same time applying a “randomization” procedure made 
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http://www.cosmosim.org/cms/simulations/bolshoip-project/bolshoip/ 
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of random translations, rotations and reflections in or der to 
avoid spurious periodicity effects dBlaizot et aljlioosl l. 

We divide a 2.4x 2.4 deg^ sky region into 200x200 pixels 
whose angular size is 43”. This corresponds to the beam size 
of a 6 m radio telescope for Ao = Acii(l + 7). The frequency 
range [238, 317] GHz {zcii =[7, 5]) is equally divided into 
60 bins with bandwidth of each bin Avo ~ 1.3 GHz. 

Given a pixel in the [G ii] map, the measured intensity 
in the frequency bin centered at i/q is 


7cii(t'o) 


1 _ ^cu 

(A0)2 47rr?(l + Zj)'^ ’ 


(19) 


where A9 is the angular size of the pixel and r-j is the co¬ 
moving distance up to Zj ; the sum is performed on all halos 
seen by this beam and with redshift 

-:-^-1 ^ ^ ^-1- (20) 

Vo + Ai/o/2 i/Q — Ar'o/2 

For the GO and [C i] emission lines the procedure is the 
same as for [C ii]. For the FIR continuum, Eq. (1191) becomes 


3.5x10^ 

„ 3.0x10^ 

" 2.5x10^ 

“ 2.0x10^ 

1.5x10^ 

1200 

^ 1000 

800 
cn 

^ 600 
400 

200 
0 

230 250 270 290 310 

Vg [GHz] 


7 0 6 5 6 0 5 5 5.0 

--- [Cll] 

- COs J 

-- [Cl(l-0)] + [ClC2-l)] 


continuum 



7FiR(r'o) = 


1 L{^SEDm(v)(l + Zj) 

( A ^)2 2 ^ 47 rr |( l -(- 2:^)2 


( 21 ) 


where v = iyo(l + Zj), and SEDih,(;z) is the normalized spec¬ 
trum template of the galaxy, 

rc/S^m 

/ SEDiR(iz)di/ = 1. (22) 

J c/lOOO^im 

For simplicity, we choose the Spi4 spir al galaxy SEP tem - 
plate from the SWIRE template librar\FI (IPolletta et al.llioOTl ) 
as a typical continuum template for all galaxies (note that 
only the FIR part is used). For all above radiation, we take 
into account the redshift distortions produced by peculiar 
motions along the radial direction. 

We also generate noise maps at each frequency by 
adopting Fays = 150 K, Apix = 128 x 128, and a total 
integration time fobs = 5000 hr. The frequency range that 
we are considering, [238 - 316] GHz, is sufficiently far from 
the prominent water atmospheric absorption at 325 GHz. 
Although there is a deep decreasing trend of transmission 
with increasing frequency, we note that the transmission is 
larger than 0.9 at vg < 300 GHz, for pwv=l mm, meaning 
that the [G ii] signal from 2 > 5.3 galaxies is not strongly 
affected by transmission issues. In case of other atmospheric 
absorption features, we could simply drop the correspond¬ 
ing frequency bins. This treatment would not strongly affect 
our conclusions. 

For the (316 ± 0.65) GHz frequency bin, maps of [G ii] 
signal, FIR continuum foreground, GO lines contamination 
and their sum plus the [C i] lines and instrumental noise are 
shown separately in Fig. [S] The GO map includes all CO 
lines from J = 1 to J = 13. Fig. [6] shows the mean [C ii] 
signal, FIR continuum, CO lines and [C i] lines as a function 
of frequency. The FIR continuum is much stronger than the 


Figure 6 . Mean intensity of [C ll] signal, CO and [C l] contam¬ 
ination and FIR continuum as a function of observed frequency. 
The ZQii is marked on upper abscissa. 


[C ii] and contamination emission: for instance, at 316 GHz 
the FIR continuum is ~ 3 x 10® Jy sr“^, while the [C ii] 
signal and the CO contamination are ^1200 Jy sr“^ and 
~800 Jy sr“^, respectively. Moreover, although at 2cii ~ 5, 
the CO lines signal is comparable to the [C ii] one, at 2 cii ~ 
7 CO dominates by a factor ~ 20. 

At zcii = 5, the sum of the two [C i] line fluxes is 
~100 Jy sr“^, therefore negligible with respect to the [C ii] 
signal. However at 2cii = 7 it represents an important con¬ 
tamination for [C ii] intensity mapping, being ~60 Jy sr“^, 
namely comparable to the [C ii] signal. 

The frequency dependence of the fluctuations of the 
[C ii], FIR, CO, ]C i] and the instrumental noise, can be 
qualitatively appreciated by inspecting a single line of sight 
cut through the mock light cone, as shown in Fig. (3 


2.6 Recovering the ]C ii] signal 


To recover the ]C ii] signal from the observed maps it is 
necessary to subtract the CO and [C i] contamination and 
FIR continuum foreground. For the purpose of contamina¬ 
tion subtraction, it is useful to know the fractional contribu¬ 
tion from halos with different mass to the total fluctuations 
signal coming from galaxy clustering. The clustering term of 
the angular power spectrum of a line emitted from sources 
located in a narrow redshift range is 


Rclust IZC 


J dMh 


(23) 


where L{M-o) is the line luminosity and b{M\y) is the halo 
- ,, , , , bias. For i>o = 316 GHz, the observed CO lines with J = 

° http://www.iasf-milaiio.iiiaf.it/-polletta/templates/swire_templatep.htinl ^ . io 

6 ,,7 -J iu- J i i n • 1 i Ji5,& and / are from 2 ~ 0.45,0.82,1.18, and 1.54, respec- 

“ We consider the third octile oi precipitable water vapour, ’ ’ > i i i r- 

pwv = 0.913 mm, namely the one assumed by the ALMA tively. Most of the CO contamination to the 316 GHz map is 

Sensitivity Calculator (ASC) in the default case. The due to these four lines. The two ]C i] lines are from 2 ~ 0.56 

reader may also refer to Fig. 2.14 of the ALMA handbook and 1.56 respectively. Fig.[8]shows the fractional power spec- 

(https://alinascience.eso.org/docuinents-aiid-tools/cycle“0/alina-tEiaimifer©lfthkai^>S)ofec^klSYdawfitotatoili^^S for [C ii] and these 
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Figure 5. Map of the [C ll] signal (left top), FIR continuum foreground (right top), the CO lines contamination (left bottom) and full 
observed map made by the sum of the signals and instrumental noise (right bottom). All maps are for the (316 ± 0.65) GHz frequency 
bin. As the emission line signal is much weaker than the continuum, the full map looks very similar to the continuum map. 


four CO lines and two [C i] lines. In order to decrease CO 
and [C i] contamination by > 90%, radiation from halos 
above ~ 10^^ —3x10^^ M© must be subtracted. In a 2.4 x 2.4 
deg^ field there are ~ 2 x 10® halos with Mh > 10^^ Mq at 
all redshifts; however, only those whose CO or [C i] emis¬ 
sion lines are received in the frequency bin centered at vo 
produce contamination. 

In principle, de-contaminating the map would be easy if 
we could measure the contamination line flux of each galaxy 
in the map, then subtract it from the relevant pixels. How¬ 
ever, this procedure is very time consuming. Alternatively, 


we could identify pixels that are supposed to be heavily 
contaminated by CO or [C i] lines and discard them com¬ 
pletely. We note that the CO or [C i] sources should be much 
brighter than [C ii] sources, say, in the optical/IR band, 
where they are also more easily resolved. Thus, one can di¬ 
rectly drop the pixels in which optical/IR bright sources are 
detected in the redshift range from which CO or [C i] lines 
are redshifted into the observed frequency bin. Of course, 
this would also cause a loss of [C ii] flux from the dropped 
pixels. However, as [C ii] (CO or [C i]) lines come primarily 
from high-z (low-a) galaxies, there are no correlations be- 
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Figure 7. The [C ii] signal, CO and [C i] lines contamination, 
instrumental noise (bottom panel), and the FIR continuum fore¬ 
ground (top panel) as a function of frequency along a randomly 
selected line of sight cut through the mock light cone. 



Figure 8. Fractional contribution from halos below to the 
clustering term (Eq. [23ll of the angular power spectrum for the 
316 ± 0.65 GHz frequency bin for various emission lines. The two 
thin lines refer to [C l(l-O)] and [C l(2-l)] lines respectively. To 
guide the eye we have plot a horizontal line corresponding to 0.1. 

tween the two galaxy populations, i.e., removing CO([C i])- 
bright pixels is equivalent to a random masking of [C ii] 
signal. The missing [C ii] p ower is limited as lon g as the 
masking fraction is ^ 30% llKashlinskv et aDl2005l l. 

How do we select optical/IR bright so urces? The natura l 
choice is to use the K-band magnitude JSilva et al.ll20lj i. 
However, this quantity may not be a good indicator for IR 
and CO or [C i] luminosities. For example, the K-band to 


total IR flux ratio of the 7 spiral galaxy templates in the 
SWIRE SED library varies by a factor of ~30. The discrepancy 
for different galaxy types is even larger. To use the specific 
K-band brightness in contamination removal, the scatter of 
the Lk — Tir relation shou ld be modele d reaso nably. 

In the samples used in iGreve et al.l (l2014h . we find 51 
(U)LIRGs and 15 high -2 DSFGs whose K-band flux can be 
found on their website. From such data, we fit the following 
Lr/ — I/iR relation: 

( „sX- ) - (tJ) + 

with a standard deviation of residuals equal to 0.35; Lr/ is 
the luminosity at the rest frame frequency that is redshifted 
into K band. For a halo with mass Mh and IR luminosity Lir 
(computed from Eq. I16II , Lr/ is randomly generated from 
a log-normal distribution with mean given by Eq. (I24II and 
standard deviation ctr/ = 0.35. We find that, for example, a 
halo with typical IR luminosity 10^^ Lq at 2 = 1 is as bright 
as niK = 21 when adopting the Spi4 spiral galaxy SED tem¬ 
plate. We are then confident that K-band luminosities can be 
safely used to remove CO and [C i] contamination. The ad¬ 
vantage with respect to relying on UV/optical magnitudes is 
that dust extinction effects in the K-band are much smaller, 
and can be neglected to a first approximation. 

The FIR foreground subtraction algorithm exploits the 
fact that the continuum is a very smooth function in fre¬ 
quency space. Such feature is widely used, for e xample, in 
HI 2 1 cm intensity mapp ing (IWang et al.l l2006l : Ijelic et akl 
l2008l: I Alonso et al1l2014b . For this reason, we believe that 
assuming the same FIR continuum template for all halos 
is acceptable, as smoothness without specifying the slope, is 
the only feature of the foreground that is required by this al¬ 
gorithm. We check that, adopting a very different SED tem¬ 
plate, e.g. an elliptical galaxy or a starburst galaxy, result in 
different slopes and amplitude of predicted FIR foreground, 
but the recovered [C ii] signal is almost identical. 

In what follows we list the steps for recovering the [G ii] 
signal from the full, observed map at the 60 frequency bins 
in which the frequency range [238, 317] GHz (2011 =[7, 5]) 
is sampled {Auq = 1.3 GHz). They form a set of 200 x 200 
line of sights (one for each pixel of the map). 

(i) Identify the “GO or [G i] contaminated” pixels (pix¬ 
els containing rriK < 22 galaxies whose contamination lines 
are redshifted into the relevant band) in each line of sight, 
replace their flux with the interpolated value from the two 
neighboring pixels along the same line of sight. 

(ii) For each line of sight, take out its foreground com¬ 
ponent that is found by either singular value decomposition 
(SVD), or polynomial fitting algorithm (details are given in 
Appendix 0). 

(iii) Set the flux of “GO or [G i] contaminated” pixels 
identified in step (i) be zero. 

After the above procedures, the final map contains the [G ii] 
signal, instrumental noise, and relatively negligible FIR con¬ 
tinuum foreground and contamination residuals. We check 
that, at angular scale ~ 1000”, for frequencies corresponding 
to 2CII = 5, 6 and 7, 0.1%, 6% and 20% of the GO contami¬ 
nation power spectrum is left as residuals, respectively. For 
[G i] lines the corresponding fraction is 2%, 5% and 10%, 
respectively. 
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8 = 2-n/q [arcsec] 

Figure 9. The remaining angular power spectrum after fore¬ 
ground and contamination removal (open squares connected with 
lines) and the original [C ll] signal (filled circles connected with 
lines) as a function of angular scale 9 = 21^jq, where q is the 
wavenumber, at zqh ~ 5 (solid), 6 (dotted) and 7 (dashed) re¬ 
spectively. We give the masked percentages for each. The assumed 
observational configuration is a 6 m telescope with Tsys = 150 K, 
^pixel = 128 X 128, and a fobs = 5000 hr. 


3 RESULTS 

Fig. [9] shows the recovered [C ii] angular power spectrum 
(dashed) from 2 ~ 5,6, 7, along with the original signal 
(solid). At 2CII ~ 5, the power spectrum is almost per¬ 
fectly recovered, with the slight deficiency at large angular 
scales due to the discarded [C ii] flux in CO or [C i] contam¬ 
inated pixels. At 2CII ~ 6, the [C ii] signal is recovered for 
9 it 1000"; at smaller scales the limiting factor is the instru¬ 
mental noise that can however be suppressed by a longer 
integration time. The signal from 2cii ~ 7 remains largely 
inaccessible, as noise dominates at all scales. 


4 CONCLUSIONS AND DISCUSSIONS 

We have studied the collective [C ii] emission signal from 
star-forming galaxies at 2 5, as well as the influence of 

FIR continuum foreground and CO and [C i] contamination 
on the experimental detection of the signal. 

To this aim we have combined the predicted [C ii] line 
luminosity as a function of the galaxy star formation rate 
and metallicity (derived from single galaxy simulations in¬ 
cluding a sub-grid treatment of the interstellar medium pre¬ 
sented in V15) with a semi-empirical approach to compute 
the Leu — Mh relation. This relation is subsequently applied 
to halo catalogs built from the large-scale N-body simulation 
BolshoiP, to generate mock maps of [C ii] signal. 

To compute the FIR continuum foreground, we derived 
the Lir — Mh relation via the abundance-matching tech¬ 
nique. As for the contamination by CO lines emitted from 
low-redshift galaxies, instead of using the poorly constrained 


CO luminosity functions, we use intermediate Leo ~ Lir re¬ 
lations that were better fitted from measurements of both 
local and high-redshift samples. We use the theoretical cal¬ 
culations of the [C i] line luminosities as a function of Lir, 
according to the Popping et al. (2014) model, mentioned 
in Sec. 2.5. We generated mock maps for FIR continuum, 
CO and [C i] emission, in close analogy with the [C ii] mock 
maps. 

We carried out FIR foreground removal and contamina¬ 
tion masking experiments on the total mock maps (contain¬ 
ing the signal -|- foreground -|- contamination and also the 
instrumental noise) to recover the angular power spectrum 
of original [C ii] maps. We pointed out that, in order to ef¬ 
ficiently subtract the CO and [C i] contamination one could 
discard pixels that are allegedly contaminated by contami¬ 
nation lines. This is feasible if the map has sufficient angu¬ 
lar resolution to avoid losing too many pixels. We estimated 
that if the intensity map has a resolution ~ 40", contami¬ 
nation can be suppressed by dropping all pixels containing 
galaxies brighter than rriK = 22 and located at the relevant 
redshift range. 

We found that the 2 > 5 [C ii] signal comes mainly 
from halos in the mass range Mq (H-band apparent 

magnitude ~ 26.8 — 23.8); as this mass range is narrow, in¬ 
tensity mapping is an ideal experiment to investigate these 
early galaxies. The [C ii] signal from 2cii ~ 5 — 6 is de¬ 
tectable for a ground-based, noise-limited telescope with a 
6 m aperture, Tsys = 150 K, a FIR camera with 128 x 128 
pixels in about 5000 hr total integration time. Although fea¬ 
sible in principle, the experiment is difficult to be performed 
using currently available telescopes. In addition, the integra¬ 
tion time itself could be longer if the atmospheric conditions 
are on average worse than assumed here. Therefore a ded¬ 
icated telescope is required and its location is essential. In 
any case, our study will serve as a robust guideline for the 
design of future facilities. 

A further motivation of a [C ii] intensity map¬ 
ping experiment is to detect the signal from faint 
galaxies that are unresolved even in deepest opti- 
cal/IR surveys. This is important as these faint galax¬ 
ies are believed to contribute m o st ionizing pho- 


20131: Choudhurv & Ferraral 2007 

: Lorenzoni et al. 2011: 

Jaacks et al. 20121: iBouwens et al. 

I 2 OI 2 I: Finkelstein et al.l 

2012 (1. 


In what follows, we therefore discuss the feasibility of 
such a kind of experiment. As can be seen from Fig.jS] at 2 ~ 
5, faint galaxies hosted in halos below 10^^ Mq only produce 
less than 1% of the total [C ii] power spectrum. Therefore 
an intensity mapping experiment aimed at detecting such 
faint galaxies could only be carried out by a telescope whose 
noise level is background-limited, such as a cryogenic space 
telescope. In this case t he noise is mainly du e to Poisson 
fluctuations of the CMB llVisbal fc Loeb|[201Cll ): 


CTJV = 


L(Tcmb, uo)hu 
X^Auqt 


(25) 


where B is the CMB emission at vq. For the following cal¬ 
culation we adopt a 2 m aperture telescope [3 and fobs = 


A 2 m aperture telescope can be considered typical for a space 
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Figure 10. The original angular power spectrum of [C ll] emis¬ 
sion from halos with [C ll] line flux < 10“^^ Wm“^ at the 
316 ± 0.65 GHz frequency bin (filled symbols) and the recovered 
one (open symbols). 

t X fimap/f^FOv = 100 hr, which is appropriate as a refer¬ 
ence for a space instrument with a low noise level. 

To analyze the [C ii] signal of these faint galaxies, the 
[C ii] flux from bright galaxies needs to be measured by a 
high-resolution interferometer array, and then subtracted in 
the relevant pixels. In the (316 ± 0.65) GHz frequency bin, 
in our light-cone there are ~ 7 x 10® halos with [C ii] line 
flux > 10“^^ Wm“^. These halos have Mh 6 x 10®® Mq. 
In addition, there are 1.3 x 10^ halos having CO and [C i] 
contamination line flux above 10~^^ Wm“^. By assuming a 
line width of 50 km s“®, to resolve these halos with a signal- 
to-noise ratio > 5, the required sensitivity is 4 x 10“® Jy. 
For comparison, at 316 GHz for a channel width of 50 km 
s“®, ALMA sensitivity with 34 antennas is ~ 4 x 10~® Jy 
with 28 hr. Assuming a ~(20 arcsec)^ FOV for ALMA, to 
cover a sky region of 2.4 x 2.4 deg^, we need ~ 2 x 10® 
pointings. Therefore, even for only one frequency bin the 
required integration time is so high (~ 5 x 10® hour!) to 
make the experiment unfeasible with current technology. 

In spite of this, we generate new maps including only the 
[C ii] signal and contamination from halos with a line flux 
below < 10“^^ Wm“^. Such halos have ^ 10®® Mq, there¬ 
fore we can use the Lcii — Mh relation in P14 simulation. We 
can recover the [C ii] power spectrum using the procedures 
detailed in Sec. mi with the following considerations. Here 
we have assumed that galaxies with [C ii]/CO/[C i] line flux 
above 10“^^ Wm“^ have already been resolved. Thus, when 
removing contamination, we can directly subtract the flux of 
the resolved sources in the relevant frequency bin, instead of 
discarding the pixel completely. Additionally, here the [C ii] 

observatory. While a larger aperture would be more helpful in 
reducing integration time, the main challenge is to measure the 
contamination flux of foreground galaxies and [C ll] flux of bright 
high-z galaxies. 


signal to FIR continuum ratio is even smaller, which allows 
us to take out the first 15 modes in continuum subtrac¬ 
tion. In Fig. [TO] we show the original and recovered power 
spectrum from zcii ~ 5. We can see that the experiment 
is indeed possible at z = 5. At higher redshift, though, the 
recovery of the [G ii] signal from faint galaxies is still ham¬ 
pered by the noise level. _ _ 

During the final stages of this work, ISilva et al.l (120141 ) 
presented a investigation of the [C ii] signal from high- 
z galaxies through mock surveys obtained from semi- 
numerical simulations. The two studies are in broad agree¬ 
ment, although they differ in the conclusions concerning the 
FIR continuum, which we found to be much stronger than 
both the [G ii] signal and contamination. Hence, accurately 
subtracting this foreground is vital in order to recover the 
[G ii] signal. The good news are that we showed here that the 
proposed algorithm to remove a spectrally smooth compo¬ 
nent from each line of sight, inspired by 21 cm experiments, 
can be successfully applied to [C ii] intensity mapping as 
well. 

Finally, we comment on a possible caveat of our work, 
related to the fact that the V15 model neglects the con¬ 
tribution of HII regions to the ]G ii] emission. We have 
already shown that HII regions in the diffuse medium do 
not have significant contribution compared with the cold 
neutral medium and warm neutral medium (see Fig. 8 in 
iGonzalez-Lopez et al.l[2014h . and VI5 further shows that 
]G ii] emission by cold neutral medium and warm neutral 
medium is negligible compared with PDRs. For what con¬ 
cerns HII regions surrounding molecular clouds (MG), we 
consider the MG properties predicted by the V15 model: 
number density, nn ~ 10® cm“®; size, tmc ~ 1 pc; ioniza¬ 
tion parameteilj U ^ 10“^®®. These fiducial valu es for wh and 
U are close to the parameters used in Fig. 2 of iNagao et al.l 
(120111) . where HII regions represent MG outskirts with col¬ 
umn density 4 x 10^® cm“^. From this figure, we can see 
that in correspondence of the MC column density expected 
by the V15 model (~ 3 x 10®® cm“®) the ]G ii] emission is 
almost 50 time stronger than the one at the edge of the HII 
regions, implying that the latter only contribute /hii ~ 2 % 
to the total emission. This calculation allows us to conclude 
that in the V15 model HII regions are expected to provide 
a negligible contribution to the total [C ii] emission, once 
compared wit h PDRs. Th is conclusion is consistent with 
the results bv lAbell 120061 ). In fact, this author finds that 
for uh = 10® cm“® and U = 10“®’®' /hii .S 10%. 

The HII regions contribution to the [G ii] emission re- 
mains, however, a very c ontroversial topic. For example, 
iOberst et all l2006l. I 2 OIII) find /mi ~ (30 - 40)% in Ca¬ 
rina Nebula, whic h is a s ingle star formation region in Milky 
Way. IVasta et al.l 120101 ) rep ort fractions spannin g a range 
5.5% — 60% in their samples. jPecarli et al.l l2014h observed 
two Lyman Alpha Emitters (LAE) at z = 4.7, concluded 
that in these LAEs (which are actually members of an in¬ 
teracting system including quasars) most of the [C ii] emis¬ 
sion is from the ionized medium. However, we also note that 


® The quoted value for U is estimated without considering the 
MC optical depth tmc to ionizing photons. If tmc is considered, 
U will become smaller and the [C ll] emission from HII regions 
will further decrease. 
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ICormier et alJ (l2015l ^ found that HIT regions contribution is 
typically less than 15% for dwarf galaxies. 

To summarize, the contribution of HII regions to the 
[C ii] emission is not clearly known, especially for galaxies 
at > 5, since [C ii] observations are still poor at these 
epochs. Given that /hii is expected to be negligible accord¬ 
ing to current state-of-the-art theoretical models, we do not 
consider the [C ii] emission from HII regions as a major com¬ 
ponent in our calculations. Accounting for this contribution 
(e.g. /hii ~ 30%) would only enhance the [C ii] power spec¬ 
trum (by a factor of ~ 2), making the signal stronger. 
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APPENDIX A: FOREGROUND SUBTRACTION 

Each line of sight is a 60-elements vector Ai where the 
smooth FIR continuum is the dominant component. Thus, 
to subtract the FIR continuum it is necessary first to find its 
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principle components, and then to subtract them from Ai. 
In other words, expressing Ai as the linear combination of 
a series of base vectors that are orthogonal to each other, i.e 
Ai = bjBj, bj indicates the contribution of each vector 
to the Ai. With descending order of bj, the summary of the 
hrst several Bj is considered as the FIR component. To find 
the principle components, the 60 x 60-elements covariance 
matrix is 



where N = 200 x 200 is the number of line of sights. By per¬ 
forming SVD on this matrix a series of eigenvectors sorted 
in descending eigenvalues Sj are got, they are Bj. Then bj 
is obtained as bj = Ai ■ Bj. Sj descends more gently af¬ 
ter a certain jmax (jmax = 2 for recovering our Fig. [O] and 
imax = 15 for recovering our Fig. HOD , therefore most FIR 
continuum is considered to be represented by the first jmax 
principle components. 

In the frequency range considered in this work, the FIR 
continuum is very close to a power law, thus we try the 
polynomial htting algorithm as well. The FIR continuum is 
found to be well fitted by polynomial up to two orders and 
the recovered [C ii] angular spectrum is quite similar. In the 
main text, we have presented results obtained by adopting 
the SVD method. However, the polynomial fitting method 
would have provided consistent results. 


